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ABSTRACT 


In  the  study  of  heating  problems  of  missiles  it  is  fre¬ 
quently  necessary  to  obtain  the  surface  integral  of  the  heat 
transfer.  This  is  usually  obtained  by  a  numerical  integration 
of  the  heat  transfer  rate  distribution.  The  purpose  of  this 
note  is  to  show  that,  for  laminar  heat  transfer  to  symmetric 
bodies  at  zero  angle  of  attack,  this  integration  can  be  done 
analytically,  if  only  a  few  reasonable  approximations  are 
made.  The  result  is  expressible  entirely  in  terms  of  quan¬ 
tities  which  also  enter  the  heat  transfer  rate  distribution  ex¬ 
pression,  so  performing  the  surface  integral  becomes  a  very 
simple  matter  of  algebraic  combination  of  known  quantities. 
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A  NOTE  ON  THE  SURFACE  INTEGRAL  OF  lAMINAR  HEAT 
FLUX  TO  SYMMETRIC  BODIES  AT  ZERO  INCIDENCE 


In  the  study  of  lieatiug  prolilonis  of  missiles  it  is  frequently  necessary 
to  obtain  llii'  surface*  integral  of  the*  he*at  transfer.  This  is  usually  oblaine*d 
by  a  numerical  integration  of  the*  he*at  transfe*r  rate  distribution.  The  pur¬ 
pose  of  this  note  is  to  show  that,  for  Lunin  ir  heat  transfer  to  symme*tric 
bodies  at  :'e*ro  angle*  of  attack,  this  inte*g  ration  can  be  done*  analytically,  if 
only  a  fe*w  re*asonable  approximations  are*  maele.  The  re*sult  is  expressible 
entirely  in  tt*rms  of  quantitie*s  which  <ilso  i*nle*r  the*  heat  transfer  rate  dis¬ 
tribution  e*xpr  e*  s  s  ion ,  so  pe*rforming  the*  surface  inte*gral  be*comes  a  very 
simple*  matte*r  of  algt*braic  combination  of  known  quantitie*s. 

The*  "local  similarity"  llu*c;ry  of  laminar  he'at  transfc'i*  to  bodies  at 
ze'ro  angle*  of  attack  gives  the*  he*at  transfer  r<ite*  as^ 


q 


r  ^  n  |x  u 
w  e 


11 


r|  w 


(T 


w 


(1) 


He*  re  r  is  the*  bocly  radius  in  the*  c  r  os  s  -  s  e*  cti  on  plane,  u^  and  11^  velocity 
and  stagnation  e*nthalpy  at  the*  e*dge*  of  the*  boundary  layer,  g  and  o\^j 

the*  de*nsity,  viscosity  and  I^randtl  numbe'r  at  the*  wall,  g-^^^y  the*  non  -  dimen - 
sienna  1  e*nthalpy  gradient  at  the  wall,  <ind  £  <i  transformed  coordinate  related 
to  the  coordinate*  x  paralle*!  to  the*  body  surface  by 
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The  inde*x  j  is  xe'ro  for  twe^- dime*  ns  iona  1  <ind  unity  for  .axis  yin  met  ric  flow. 


T'o  inte*g  r.i  t  e*  c\  ove*r  the  boely  surface*  up  to  st.ation  we  multiply  by 
the*  e*le^me*nt  of  surface  .i  re*a  rlA  =  4^(7:  r)J  dx  anel  inte'grate*  from  x  =  0  to 
X  -  x^,  obtaining 


CiclA  =  Z(tt  )j 
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It  is  this  inte'gral  which  we*  will  .appre)ximale'. 


First  noti?  that  11^,  is  a  cons ta at.  S(*toa(l,  note  from  Fc|.  (^)  that  the 
numerator  of  the  t(*rm  in  tiu*  scjuare  brackc*ts  is  fl  $  / clx  so  wc?  can  write 

Q;  =  i  (,  )j  II,.  r*  ^ 

•/  o  O’ 

w 

The  Francltl  number  may  ho  t<iken  to  b(*  a  constant,  o.  Tlu*  non-dimensional 
enthalpy  giMdient  floes  not  vary  much  <iround  tlie  l)o<Iy  (‘ither^,  since  it 

has  a  weak  dependence  oti  piu^ssiirc^  eradifuit  and  a  roughly  scjinirc*  root  dtj- 
pendiMice  on  pp  ratio.  In  fact,  Lees^  took  it  to  be  a  constant  at  tin*  Z(‘ro 
pressure  gr.idient,  (  cuistant  ,  value  in  his  lurat  transf(‘r  tht‘ory.  Thus,  a 

re.isonable  approximation  is  to  take*  <i  suitalde  average  value,  say  g^^vvo* 
may  then  integrate  Kc\,  (-1)  to  oI)tain 

g  J  li  (xj  /  (T.  (5) 
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This  formula  (‘xpress(‘s  th('  surface  ii-itt^gral  of  tlu^  havil  transf(‘r  rate  in 
terms  of  the  transformed  bc.tdv  coordinaU*  ^  at  the  end  station  of  the  integra¬ 
tion  interval  xp  lIoweV(*r,  ^  can  l>e  ('\p  massed  in  t(‘rms  of  cj  fror^i  F.cj.  (1), 
SO  cai;  be  r (da ted  to  bv 

Q,  =  ^  (  TT  r,)-'  ir  / 


\  more  iKscdul  lorir'.  lor  blunt  borlies  is  oOtiinod  In'  introducing  the 
stagnation  point  le,it  iMte  !)v  th(‘  i  pp  in)|)  r  ia  t  (‘  limiting  lorm  ol  (1): 
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'.’■lu*r(‘  (flu^,  ax)  i.‘  the  .'tagnition  pi;int  \  <do(itv  graflient.  I  nde  r  t  h(‘  a  pj:)  r  ox - 
i'nation  tliat  constant,  divisicui  id  1>:.  (n)  bv  I'.cj.  (7),  and  use  of 

:  g .  (  I  )  for  (\  \- 1  e  1  d  s 
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wl.i'.i  Ai  1.  tie  body  i  loss  .oction.il  o  i*ee  ot  x  \p  I'ld  tin*  princ  ipal 

rf‘*  Uii  ol  the  present  note,  -I'.owing  that  tlu‘  surf, ice  intt‘ginil  oi  cj  is  in- 

cer.'t'i,  f)  ropo  rt  ioiia  1  to  tin*  \.  due  ol  at  th(‘  imuI  point  ot  tlie  integr.ition 

interc.d.  Since  tlie  total  Iieat  input  to  tlie  body  up  to  station  x^  is  p  r  oport  ion.i  1 
to  tilt  )  indar\  I.iyer  thukness  tlu^re,  v dm  h  in  turn  is  itu'ersely  proportional 
to  tlie  b  i  .  I  lun  t  traiishM*  I’atc,  the  result  iii  Im].  (8)  is  phvsic  ally  n  (‘,i  s  oinible . 
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Kurth<*r  simplification  of  Kc\,  (H)  is  obtained  if  the  Imdy  can  be  taken 
to  liavc'  a  ('ool  wall  at  which  tlie  p(*rfect  p;fis  law  is  apjil  ica  bh* .  Then  we  have* 
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wlierc*  p  is  the  pressure  in  the  boundary  layer.  Finally,  if  the  wall  can  b(* 
t.ikenas  constant  te  in  j)e  r.i  t  u  r  c* ,  the  factv)r  in  parenth(*sis  in  Fq.  (9)  is  unity 
and  we  have  the  simple  resuP 
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As  rin  example*,  for  a  he*  n"!  i  s  phe  r  ic  a  1  nose  of  radius  at  hypersonic 
speeds,  iq.  is  practically  a  linear  function  of  x  =  HP  ,  an  !  o  r  =  R  sin  0  . 
The*  are*a  <ind  \'e*locity  f.ie  tor  is  then 


A  .  u  .  r .  . 

^  ^  =  TT  ir  0  sin  0  . 
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In  the*  con.stant,  low  t(*mperature  wall  case,  Ref.  Z,  Fcj.  (15)  or  Ref.  3, 
K<.\,  (7)  (with  the  correction  cjf  <i  sepiare  root  ii'i  tin*  de noniinat o r )  gives 


0  sin  0 


''<-s  [K,  (0)  M^ia  (9)] 


l/Z 


1  -  I  ^  sin""  P  , 


(  l^) 


F  (P)  0"^  -  0  sin  ^  sin^  P, 


F,  (0)  -  i/l 


i)^  -  0  sin  li;  (  1  +  Z  (sin""  0)/ ^  sin""P(I  t  (sin""P)/3) 
Thus,  the  integratc'd  heat  transI(M*  become* s,  from  Paps.  (lO)-(lZ), 
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a  result  wliich  agree*s  with  a  dire(  t  integration  of  Paj.  (IZ). 


If  a  direct  integral  iorinulation  is  desired,  instead  of  an  t*x|:> ress ion 
in  terms  of  q,  we  can  restore  f  in  place  of  q  in  K(\.  (8)  by  means  of  Ecjs. 
(1)  and  (7): 
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Further  simplification  for  constant,  low  temperature  wall  leads  to  the 
formula  analogous  to  Fa\.  (10): 
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If  the  hemispherical  nose  example  of  Eqs,  (11)  and  (  1  Z)  is  used  in 
Eq.  (15),  the  integration  can  oe  performed,  and  the  result  agrees,  of 
course,  with  P'cj.  (13). 

The  atlvantagt^  of  Eqs.  (8)  and  (10),  over  Eqs.  (1*1)  and  (15),  is  that 
the  integrals  in  the  latter  art*  dir<.M,:tly  ri*lated  to  the  heat  transft*r  rate  dis¬ 
tribution,  and  so  have  already  bi*tMi  calcuhited  if  tliat  distribution  is  known. 
In  this  case  tlu*  surTica*  intt‘gral  can  be  found  from  Ec|s.  (8)  and  (10)  with  no 
further  integration.  If,  hov\a'Ver,  on»'  flesirt's  to  find  the  integ  ratted  heat  in¬ 
put  without  first  finding  the'  heat  flux  distribution,  l-a}s.  (!•!)  and  (15)  are  the 
us(.‘ful  ont‘S. 

If  the  integral  of  q  itself  is  ch*  sired,  the  va  liu*  of  should  be  found 
from  the  formulas  given  in  Kef.  *1. 
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